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In this article, we calculate the form factors and the coupling constant of the 
gB*^B*K vertex in the framework of the three-point QCD sum rules. Three point 
correlation functions responsible for the vertex are evaluated by considering both B* 
and K mesons as off-shell states. The form factors obtained are different if the B* 
or the K meson is off-shell but give the same coupling constant. 
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I. INTRODUCTION 

During the last decade, both theoretical and experimental studies on heavy mesons 
have received considerable attention. Studying the properties of bottom-strange orbitally 
excited states can provide tests for various models of quark bound states[l|. Recently, two 
orbitally excited narrow Bg mesons are reported by the CDF collaboration with masses 
mB,^ = 5829.4 ± 0.7 MeV and m^.^ = 5839.6 ± 0.7 MeV fl, which are assigned as the 
= (l"*", 2"*") doublet states in the heavy quark effective theory The DO collaboration 
confirmed the B*2 state with mass m^*^ = 5839.6 ± l.l(stat.) ± 0.7(syst.) MeV in fully 
reconstructed decays to B^K~ In contrast, the Bg states with spin-parity = 
(0"*", 1"*") still lack experimental evidence because they are supposed to be broad decaying 
through an iS-wave transition. 

Based on the potential quark models, heavy quark effective theory, and lattice QCD 



10|, masses of the Bg mesons with (O"*", l"*") have been estimated. With QCD sum rules 
Wang studies masses of the bottomed (0+, 1'^) mesons. The calculation indicates that 
the central values are below the corresponding BK and B*K thresholds respectively, which 
is just like the case of their charmed cousins DgQ and D*^. This prediction is compatible 
with speculation in Refs. j6|, |l2|. For this reason, the strong decay of B*-^ — )• B*K is not 
predicted to occur. As a consequence, the width of S*]^ is very narrow since the main decay 



channel for this state is isospin symmetry violating B*^ — )■ B*ri — )■ B*it^. In Ref. 13|, the 
isospin-violating strong decay width of B*^ is calculated, which is around several tens of 
keV. 

Knowledge of the heavy-heavy-light mesons vertices is of crucial importance to esti- 
mate the strength of the hadron interactions. They are fundamental objects of low energy 
QCD. On the experimental side, the excited Bg mesons will be copiously accumulated at 
the CERN LHC, which makes investigations of their decay characters promising. Their 
coupling constants can bring important constraints in constructing the meson-meson poten- 
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tials and help us to analyze results of existing experiments in the framework of the meson 
exchange model 1^. For example, Faessler et.al used a phenomenological Lagrangian ap- 
proach to calculate the strong and radiative decay of the Bg doublets with = (0"*", 1"*"), 
which were considered as bound states of the BK and B*K mesons respectively 15|. In 
the process of their calculation, the coupling constant of B*^B*K vertex was an necessary 
input parameter. However, such low-energy hadron interaction lie in a region which is very 
far away from the perturbative regime, precluding us to use the perturbative approach with 
the fundamental QCD Lagrangian. Therefore, we need some non-perturbative approaches. 



such as QCD sum rules [16|, ll7| , to deal with processes relevant to hadron physics. Taking 
the point of view that the bottomed (0+, 1+) mesons Bgo and B*^ are the conventional bs 
mesons, strong coupling constants QBsoBK and gB*_^B*K are calculated with the light-cone 
QCD sum rules 18|. The bottomed 0^ meson BgQ is identified with the conventional bs 



19|. 



meson and Bg^BK vertex has been calculated in Ref. 

In this article, the form factor and the coupling constant of the gB*-^B*K vertex is cal- 
culated in the framework of the three-point QCD sum rules. In the calculations we use 
the same technique developed in the previous works for the evaluation of the coupling 
constants in the vetices D*Dtt HQ) DDp^, D*D*tt DsD*K, D*DK [24], DDu 
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25|, D*D*p (26|, DgjDK [27| and D^DK* 

This paper is organized as follows. In Sec. Ull we give the details of the QCD sum rules 
for the considered vertex when both B* and K mesons are off-shell. Sec. Illll is devoted to 
the numerical analysis and we compare them with results obtained in other works. 



II. THE SUM RULE FOR THE B*^B*K VERTEX 

In this section, we obtain QCD sum rules for the form factor of the B*^B*K vertex. 
The three-point function associated with the B*^B*K vertex, for an off-shell B* meson, is 
given by 

r?:jp,y) = I d'xdS e^^'-e-^-^(0|T{jf (x)jf (2/)j,^-t(o)|o), (1) 

where the interpolating currents are j^ix) = s{x)ja'l5u{x), j^* (x) = u{x)ji,b{x), and 
jfj-'^ix) = s(x)7^ 756(0;). The correlation function for an off-shell K meson is 

rf.b,p') = I d'xd'y eV-e-.-. (o|T{jf (x)/'(2/)j^^-t(o)}|o) , (2) 

with the interpolating currents j^' [x) = u{x)'jiyb{x) , j^{x) = s{x)'y5u{x), and jl^^^{x) = 
s(x)7^ 756(0;). In the expressions, q = p' — p is the transferred momentum. 

The general expression for Eq. ([1]) has fourteen independent Lorentz structures, which 
can be decomposed in terms of the invariant amplitudes associated with each of these 
tensor structures in the following form: 

^^Lip^p') = ^l(^)^^)'^ q^)gf,uPa + ^2(p^p'^ q^)g^,aPv + r3(/,p'^ q^)guaPt, 
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+r7(p^p'^ q'^)p^,PuPa + Tsip'^^p''^, q^)p'^PvPa + Tq{p'^,p''^, q'^)p^,p'^Pa 
+^lo(p^p'^ q'^)Pf,p^p'a + q'^)p'f,plPa + ^l2(p^p'^ q'^)p'f,PuP'a 

+^l3(p^p'^ q^)Pt,p'uP'a + ^l4(p^p'^ (t)p>i,p>vP>c.- (3) 

In principle, we can work with any Dirac structure appearing in Eqs. Q. However, there 
are still some points which one must follow: (i) The chosen structure must appear in the 
phenomenological side, (ii) The chosen structure should exhibit good OPE convergence, 
(iii) The chosen structure should have a stability that guarantees a good match between 
the two sides of the sum rule. After some calculations, we find that the structures which 
obey these points are p'^p'yP'a in the case of B* off-shell, and p^^J^p'y in the case of K off-shell. 

In order to get the sum rules, the correlation functions need to be calculated in two dif- 
ferent ways: in phenomenological or physical side, they are presented in terms of hadronic 
parameters such as masses, leptonic decay constants and the form factor; in theoretical or 
QCD side, they are evaluated in terms of QCD degrees of freedom via operator product 
expansion (OPE) in deep Euclidean region. The sum rules for the form factors are ob- 
tained with both representations being matched, invoking the quark-hadron global duality 
and equating the coefficient of a sufficient structure from both sides of the same correla- 
tion functions. To improve the matching between the two representations, double Borel 
transformation with respect to the variables, = —p^ — and P'^ = — p'^ — ?■ M'^, is 
performed. 

In detail, we calculate the physical part of the first correlation function ([T]) for an off-shell 
B* meson. The form factor gB*-^B*K{q'^) is defined by the matrix element 

{B:,{p)\K{p')B*iq)) = -igB:,B'K{q')ri* ■ e, (4) 

where e and r] are the polarization vectors associated with the B* and 5*^, respectively. 
The meson decay constants /b*^ , /b* , and fx are defined by the following matrix elements: 

(0|tf-|i?:i(p)) = ^biJbi^V,. 

(0|jf = ifKPc (5) 

Saturating Eq. ([1]) with B^^, B* and K states, using Eqs. (jl]) and ([5]), and then summing 
over polarization vectors via 

* - _ I Mil 

Vf^Vu = -9f.u + (6) 
the physical side of the correlation function for B* off-shell can be obtained as 
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m 



A similar procedure is carried out to obtain the final expression of the physical side of 
the correlation function for an off-shell K meson as 



sl 



in which represents the contributions of the higher states and continuum. 

In the following, we turn to the QCD side of the correlation functions in deep Euclidean 
space. For an off-shell B* meson, each invariant amplitude T^{p' ,p) appearing in Eq. ([3]) 
can be written in terms of perturbative and non-perturbative parts: 

= (n;,, + n^„^„^,,) . (9) 

In the expression, the perturbative part is defined as 



where pi{s, u, Q"^) is called spectral density. The spectral density is obtained by calculating 
the bare loop diagrams (a) and (b) in Fig. ([T]) for B* and K off-shell, respectively. We 
calculate these diagrams in terms of the usual Feynman integral with Cutkosky rules, 
i.e., by replacing the quark propagators with Dirac delta function ^21^2 {—2-ni)5{q^ — 
m?)6{q^). After some straightforward calculations, the spectral density associated with the 
structure p'^p'^p^ is obtained as following: 

p^* (s, t, u) =--^^— _[m^(s=' - 3s2(f - 3m) + ?>s{e - Atu + Bm^) - (t _ uf) 

+ m^(3m^(3s^ + s^{'iu - 5t) + s{t^ + Atu - hu^) + (t - uf) - 2/ + s^(5t - 13m) 

+ s^{-?,t^ + Atu + 3^2) - s{t - llu){t - uf + {t- m)^) 

+ m^(3m^(3s3 + s^{M - hu) + s{-ht^ + Atu + u^) - {t - u)^) 

- 2ml{As^ - s^{t + u) + s^{-9t^ + 22tu - 9^^) + 5s{t - u)^{t + u) + {t - u)^) 
+ s(/ - s^{t - 5u) + s^{~3t^ + lOtu - 9u^) + s{5t^ - Ut^u + Utu^ - u^) 

- 2{t - u)^{t + 2m))) + ml{s^ + s^{9t - 3m) + 35(3^^ - 4tM + u^) + {t - m)^) 

+ m^(-2/ + s^(5m - 13t) + 5^(3^2 + 4tM - 3u^) + s{t - u)^{llt - u) + {t - u)^) 
+ mls{s^ + s%5t -u) + s'^{-9t^ + lOtu - 3u'^) - s{t^ - Ut^u + lltu^ - hu^) 
+ 2{t - M)^(2t + m)) + s^{s^{-{t + u)) + s^{3t^ - 2tu + 3u^) 
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+ s(-3r + t'u + tu' - 3m^) + (t - uY{r + Atu + u^))] (11) 
With a similar procedure, we obtain the spectral density 

t, «) = [X(7^77jpK(^' - 2^(5 + ^) + + + 

+ m^(m^(4s2 - 2s(f + u) - 2{t - uf) - {s - t + u){s'^ - 2u{s + t) 
+ Ast + + u^)) + mtii-s - tf + Asu - 2tu + u^) 
- ml{s + t-u){{s- tf + Asu - 2tu + u^) 

+ s{s^{t + m) - 2s{f -tu + M^) + (t - M)^(t + u))] (12) 

for the off-shell K meson associated with the structure p^p^. All powers of the strange 
quark mass are included when calculating the spectral density. Herein A(a, b, c) = + 6^ + 
- 2ac - 2bc - 2ab and t = = -Q^. 

The nonperturbative contributions in the QCD side are determined from the quark 
condensate diagrams (c), (d), (e) and (f) of Fig.([T]). As what has been shown in Refs. 
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20|] , heavy quark condensate and gluon condensate contributions are negligible as compared 



with the perturbative one. As a result, only light quark condensates contribute to the 
calculation. It is noticed that contributions of diagrams (d), (e) and (f) are zero after the 
double Borel transformation with respect to the both variables and P'^ . Hence, we 
calculate the diagram (c) for the off-shell B* meson and obtain 

for the off- shell B* meson. 

We use the quark-hadron duality assumption to subtract the contributions of the higher 
states and continuum, i.e., it is assumed that 

phigherstates^^^ u) = p°^^(s, U, t)6'(s - SQ)e{u - Uq), (14) 

where sq and uq are the continuum thresholds. 

The next step of the method is to apply the double Borel transformation with respect 
to the P2 = -p^ and P'^ = -p'^ M'^ to the physical as well as the QCD side. 

We get the final sum rules for the corresponding form factors as 
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A 7r2 



ds I dup {s,u,t) 



(15) 



and 



„2 



K (n^. 'imsmB*{Q^ + ml) ^ 
9b'B'k{Q = -f — T-f 2 e^e^ 



-1 
4^ 



ds I dup^ {s,u,t) 
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B* 



K 




Bti (e) K B;, (f) B* 

FIG. 1: (a) and (b): Bare loop diagram for the B* and K off-shell, respectively; (c) and 
(e): Diagrams corresponding to quark condensate for the B* off-shell; (d) and (f): Diagrams 
corresponding to quark condensate for the K off-shell. 



CM^eM^ (16) 

for the off-shell B* and K meson associated with the B*^B*K vertex, respectively. In 
the above expressions, we have Umin = itT'I ~ ^ind Umax = ^ + t ~ ~ 
integration regions are determined from the fact that arguments of the three 5 functions 
must vanish simultaneously with consideration of the three Heaviside step functions. 
It is noticed that in this work we use the relations between the Borel masses and 

2 

with ^ = ^4f^ for a B* off-shell and ^ = ^ for a K off-shell. 

M' 0.64 M'^ mt,f 
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M2[Gev2] (a) 




FIG. 2: a) The dependence of the form factor g^l B'xiQ^ ~ GeV^) on Borel mass parameters 
for Ag = 0.5 GeV and = 0.5 GeV and b) pole-continuum contributions. 



III. NUMERICAL ANALYSIS 



This section is devoted to the numerical analysis of the sum rules for the form factor. 
Input parameters are shown in Table [B The values /s*^ , t^b*-^ and are used the same 
as that in Ref. We first determine the three auxiliary parameters in the sum rules, 
namely the Borel mass parameter and the continuum thresholds, sq and uq. The 
continuum thresholds, sq and uq, are not completely arbitrary as they are correlated to the 
energy of the first excited states with the same quantum numbers as the state we concern. 
They are given by sq = ("^s'j + A^)^ and uq = (m + A„)^, where m is the K meson mass 
for the case that B* is off-shell and the B* meson mass for that K is off-shell. Since these 
parameters are not physical quantities, we need to look for the working regions at which 
pole dominance and stability of the sum rule with the Borel mass parameter are satisfied. 
Using = Au = 0.5 GeV for the continuum thresholds and fixing = I GeV^, we find 
a good stability of the sum rule for gB\B*K with > 10 GeV^, as is shown in Fig. (Et). 
Fig. ^p) demonstrates that contributions from pole terms with variation of the Borel 
parameter M^. We see that the pole contribution is larger than continuum contribution 
for M2 < 17GeVl We choose = ISGeV^ as a reference point. 



TABLE I: Parameters used in the calculation. 



m,(GeV) mb(GeV) mB*^(GeV) mB*(GeV) mx(GeV) /B*j(GeV) /B*(GeV) /i^(GeV) (ss)(GeV)' 



0.13 



4.7 



5.72 



5.325 



0.49 



0.24 



0.17 



0.16 0.8 * (-0.23)2 



Now, with the determined and other input parameters, the dependence of the form 
factors on is plotted in Fig. |3l The circles correspond to the g^l b*jcLQ11, form factor 

si r I 



in the interval where the sum rule is valid. Following previous works 19|, |22|, we fit the 
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FIG. 3: g% (circles) QCDSR form factors as a function of Q^. The solid line correspond 
to the monopolar parametrization of the QCDSR data. 

circles with B* off-shell by the monopolar parametrization, 

_ 1739 GeV3 

The coupling constant is defined as the value of the form factor at = —m^, where m is 
the mass of the off-shell meson. Using = — m^, in Eq. (fT7|) . the coupling constant is 
obtained as gB\B*K — 24.96 GeV. 

In the case that K is off-shell, Fig. (jUi) demonstrates a good stability of gB*^B*K with 
respect to the variations of Borel parameters for > 15GeV^. We see that the pole 
contribution is bigger than the continuum one in the Borel window < 25 GeV^ from 
Fig. (jib)- We choose = 23 GeV^ as a reference point. Our numerical results can be 
fitted by the exponential parametrization 

9Bt,B'Km = 19.48 Exp (p^^l^) GeV, (18) 

shown by the solid line in Fig. O Also, g^* = 22.29 GeV is obtained at = —Tn\ in 
Eq. f lTSj) . which is in a good agreement with the result g^l o.^ = 24.96 GeV. Taking the 

si 
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FIG. 4: a) The dependence of the form factor g^* ^t^iQ"^ — 1-0 GeV^) on Borel mass parameters 
for As = 0.5 GeV and A^, = 0.5 GeV and b) pole-continuum contributions. 



35 1 1 1 1 1 1 1 1 r 




Q^[GeV^] 

FIG. 5: c/f. B*KiQ^) (circles) QCDSR form factors as a function of Q^. The solid line corre- 

si 

spond to the exponential parametrization of the QCDSR data. 

average of these results, we get 

gB*^B*K = (23.62 ± 1.33) GeV. (19) 
Comparing Fig. [3] with Fig. [Sj we also discover that the heavier is the off-shell meson, the 
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□■^[GeV^] (a) Q^IGeV'^] (b) 



FIG. 6: Dependence of the form factor gBi^B'KiQ"^) on the continuum threshold for both 
B* and K off-shell cases. The dotted curve corresponds to A^^n = 0.4GeV, the solid one to 
As,u = 0.5GeV and the dashed curve to /S.s,u = 0.6 GeV. 

TABLE II: Theoretical estimations of the strong coupling constants from different models 

Coupling constant This work [30, ^Ij [181 
gBi^B'KiGeV) 24.53 ± 5.20 23.572 18.1 ±6.1 



more stable is its form factor in terms of Q^, which was first observed in Ref. 22 1. 

In order to study the dependence of our results on the continuum threshold, we vary 
A^ u between 0.4 GeV < A^ ^ < 0.6 GeV in the parametrization corresponding to the case 
of an off-shell B* and K. As can be seen in Fig. |6l this procedure gives us an uncertainty 
interval of 20.55 GeV < q^I^b^k < 29.73 GeV and 19.33 GeV < g'^*^B*K < 27.87 GeV. 
Taking into account uncertainties due to the continuum threshold parameter, we finally 
obtain 



gBi^B*K = (24.53 ± 5.20) GeV. 



(20) 



Together with the predictions from other theoretical approaches, the numerical results of 
the coupling constant are presented in Table [Tll Based on the heavy-light chiral lagrangian, 
the unitarized two-meson scattering amplitudes are analyzed, and the strong coupling 



constants gB*-^B*K = 23.572 GeV is obtained [30|, l3l|- Our result is in good agreement with 



their estimate. Whereas, it is estimated to be 18.1 ± 6.1 GeV with light-cone QCD sum 
rules 18|, the central value of which is about 30% smaller than above two predictions. 

In conclusion, we have used three-point QCD sum rules to calculate the form factor of 
the B*iB*K vertex. Both cases that B* is off-shell and K is off-shell have been considered. 
As a side product of the form factor, the coupling constant is estimated, which is compatible 



with the results in Refs. 30|, |3l| but larger than that from light-cone QCD sum rules [l8|. 
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